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Abstract. On the basis of Bethe ansatz sohition of two-component bosons 
with SU(2) symmetry and (5-function interaction in one dimension, we study the 
thermodynamics of the system at finite temperature by using the strategy of 
thermodynamic Bethe ansatz (TBA). It is shown that the ground state is an isospin 
"ferromagnetic" state by the method of TBA, and at high temperature the magnetic 
property is dominated by Curie's law. We obtain the exact result of specific heat and 
entropy in strong coupling limit which scales like T at low temperature. While in weak 
coupling limit, it is found there is still no Bose-Einstein Condensation (BEC) in such 
ID system. 
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A two-component Bose gas has been produced in magnetically trapped ^^Rb by 
rotating the two hyperfine states into each other with the help of slightly detuned Rabi 
oscillation field It was noticed that the properties of such Bose system can 
be different from the traditional scalar Bose system once it acquires internal degree of 
freedom. Bethe ansatz solution of SU(2) two-component bosons in one dimension was 
obtained]^, It was pointed out that the ground state of such a system is an isospin 
"ferromagnetic" state which differs from that of spin- 1/2 fermions in one dimension, 
such as SU(2) Hubbard model p etc.. 

An interacting SU(2) boson field trapped in a one dimensional ring of length L can 
be modeled by the following Hamiltonian 



H = dx 



(1) 



where a,b = 1, 2 denotes the z-component of isospin. The Bethe ansatz equations (BAE) 
of eq. are obtained as follows 



kj — ki + ic kj — \u — ic/2 
, ki - ki-ic^\ k^ - X^ + ic/2 

1 = 1 1^=1 



^ _ _ -TT Xy- ki- ic/2 -pr - + ic . 

il + %cl2 \-\,-ic ^' 

where M denotes the total number of down isospins. Eq. (|^) differs from the BAE 
of scalar bosons with periodic conditional . The second equation of eqs. (0) of isospin 
rapidity A arises from the application of quantum inverse method]^, which can be 
inferred from spin-1/2 fermionsQ] too. However, the symmetry of bosons wave function 
leaves the first term on the right side of first equation of eqs. (0), which does not appear 
in the BAE of fermions. 

The strategy we use here is the thermodynamic Bethe ansatz (TBA) which was 



pioneered by C. N. Yang and C. P. Yang for the case of the delta-function Bose gaspO 
It is used to derive a set of nonlinear integral equations called TBA equations, which 
describe the thermodjTiamics of the model at finite temperature. Moreover, the As can 



be complex roots which should form a "bound state" with other As|TT|] when T 7^ 0, 
which arises from the consistency of both sides of the BAE. For ideal A strings of length 
m the rapidities are A"-' = A^ + + 1 — 2f)iu + 0(exp(— 5A^)). Here u = c/2, a 
enumerates the strings of the same length m, and j = 1, . . . , m counts the As involved 
in the ath A string of the length m, A" is the real part of the string. 

Taking logarithm of the BAE (H) by using string hypothesis we arrive at the 
following discrete Bethe ansatz equations 

27CI, = k,L + 2 J2 02(% -ki)-J2 ^n{k, - a:) 

I an 

2nJ: = 2 J2 0n(A: - A:,) - 2 ^ A^uQtiX: - Xi) (3) 

I bl,t=iO 
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where 0„(x) = tan ^{x/nu) and 

{1, for t = n + l^\n — l\ 
2, fort = n + /-2,-- -,1^-/1 +2 
0, otherwise. 

Ij and play the role of quantum numbers for charge rapidity and isospin rapidity 
respectively. In order to guarantee linearly independent of wave function, all quantum 
number within a given set of {/} as well as that in {J} should be different. An 
arbitrary quantum number may be in a set (charge roots) or missing (charge hole). 
In thermodynamic limit, the distribution of charge rapidities becomes dense and it is 
useful to introduce the density function for charge roots and holes respectively. We 
denote with p{k) and p^{k) the density function of charge roots and holes, in a similar 
way, with o"„(A) and cr^ the density function of n-strings roots and holes on real axis. 
They are defined by 

p{k)+p''{k)= {l/L)dI{k)/dk 

aniX) + a^(A) = il/L)dr{X)/dX. (4) 
Then from eqs. (P) we obtain a set of coupled integral equations. 



P + p'= ^ + K2{k)*p{k)-Y,Kn{k)*an{k) 



271 

n 

a„ + = ir„(A) * p(A) - A^uKtiX) * ai{X) (5) 

where Kn{x) = nu/'K{n'^u^ + x^), and * denotes a convolution. 

In terms of the density functions of charge and isospin roots, the kinetic energy 
per length has the form E^/L = J k'^p{k)dk, the total number of down isospins is 
M/L = J2n ^ I '^n(A)(iA and the particle density of the model is D = N/L = J p{k)dk. 
If we consider the energy arising from the external field VL which is the Rabi field in 
two-component BEC experiments, the internal energy of the model is 

E/L= 1 {k^ -n)p{k)dk + ^2nn 1 GndX. (6) 

n 

And with the help of the approach first introduced by Yang and Yang[|To[|, the entropy 
of the present model at finite temperature is 

S/L = J[{p + p^) ln(p + p'^) - p In p - p'^ In p^fc (7) 

The Gibbs free energy of the model then is defined hy F = E — TS — pN, where p 
is the chemical potential. In order to obtain the thermal equilibrium, we minimize the 
free energy with respect to all the density functions subjects to the constraint (H). In 
addition, the total number of particles, the magnetization are to be keep constant. For 
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this purpose, the chemical potential /i and external field Q play the role of Lagrange 
multipliers. 

It is useful to define the energy potential for charge sector and isospin sector: 

K{k) = e^'^^l^ = p''{k)/p{k) 

r/„(A) = a^(A)K(A). (8) 

Applying the minimum condition 6F = gives rise to a revised version of Gaudin- 
Takahashi equations 

Tlnn = e{k) = k'^ - p-VL- TK2{k) * ln[l + k'^] 
-TY,Kn{k)*\n[l + r^-^] 

n 

Inr/i = ^sech(7rA/2M) * ln[(l + k'^){1 + r/2)] 

ln?7„= ^sech(7rA/2M) * ln[(l + ?7„_i)(l + r7„+i)]. (9) 
And these equations are completed by the asymptotic conditions 

lim [ln?7„/n] = 2x (10) 

n— »oo 

where x = Q/T. Eqs. (^) can be solved by iteration. Note that eqs. (j^) together with 
eqs. completely determine the densities of charge roots and isospin roots in the state 
of thermal equilibrium. The Helmholtz free energy F = E — TS is given by 

T/ r 

F = - — j + e-''^]dk. (11) 

The above approach called TBA is universal for discussing the thermodynamics 
of one dimensional integrable model. Once eqs. (|^) are solved, all thermodynamic 
quantities can be obtained from eq. (^) in principle. 

Magnetic property: zero and high temperature limit: The state at zero 
temperature is the ground state. The Fermi surface is determined by e{kp) = 0. Since 
there is no hole under Fermi surface, we can take the energy potential k = p^/ p as zero. 
As a result, from eqs. (|), it is easy to see rjn 00, and M = 0, the "ferromagnetic" 
ground state. The first equation of eqs. (P) becomes 

eo{k) = e - p - n + K2{k) * eo{k) (12) 

which gives the solution of dressed energy, and the ground-state energy may be given 
in terms of eo 



1 



kp 



Eo/L = — / eoik)dk. (13) 

Consequently, the ground state of ID SU(2) bosons is an isospin "ferromagnetic" 
state, which coincides with the analysis of Li et al.[0. Then the property of the model 
at T = is the same as that of scalar bosons in one dimension which has been discussed 
extensively by Lieb and Liniger^l. In the isospin space, however, the SU(2) symmetry 
of whole system around the ground state is broken. 
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In the high temperature hmit T — oo (free isospins), however we can assume that 
all functions ?7„(A) are independent of A. Then eqs. (0) can be written as follows, 



'1 + ?7„_i)(l + r7„+i) 
where we have neglected the term (1 + 



in the second equation of eqs. 



solution of rjn are then constants fixed by the field boundary condition (^) to be 



Vn 



sinh(n + l)a; 
sinh X 



- 1. 



(14) 
The 

(15) 



After perform the Fourier transformation on eqs. (^), we get the solution of the 
densities of A n-strings, 



— sech[7rA/2'u] * [p + a 



an + a. 



4:U 



sech[7rA/2n] * [a. 



h 

n+l 



(16) 



If we assume that (j„ and are independent of A or let c = 0, the total number of down 
isospins has the form, 



E 



P 
2 



-CTr, 



(17) 



where is maximal length of A strings. In the absence of Rabi field, we have 
M/N = 1/2, the system at high temperature is a quasi "paramagnetic" state. If the 
external field Q is small, expanding eq. (^) for small filed x and integrating the equation 
over A space, we get the magnetization of the model, let be the total number of 
isospin rapidities in all n^-strings. 



L 



Mm 

2L 



+ 



:i^ 



T ■ 2T2 

where the first term in the brackets arises from self-magnetization, while the others are 
contributed by Rabi field. Eq. (0) indicate that the magnetic property of the model 
in high temperature regime dominated by Curie's law xool/T. 

Strong coupling limit: When r/ — > oo, Kn{k) = 0, from eqs. (^ we have 

e = e-n-fi. (19) 

The free energy of the system ( pA]) at low temperature now can be solved by step 
integration. 



F/L = fiD 



3^ ^ 24/xV2 



(20) 



where the external field is set to zero. 

We can not deduce the specific heat directly from the free energy obtained above 
because the chemical potential is a function of temperature. From eqs. (^, the density 
of charge rapidity has the form 
_ 1 1 

^~ 2^1 + e(^'-A')/^" ^^^^ 
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Clearly, at zero temperature, the Fermi surface is just the square root of the chemical 
potential, so we have fiQ = tt'^D^. At low temperature, however, it is determined by 
D = N/L = J p{k)dk. After integration, we have a time dependent chemical potential 



1 



(22) 



Then the free energy becomes 



F/L = fioD 



1 + 



3tt 



3/2 



TdS/dT, we find the specific heat 



1 + 



(23) 



Since by thermodynamics S = —dF/dT and 
at low temperature is Fermi-liquid like 



S 



T 



(24) 



It is the same as the result of one-component case, since for the strong coupling limit 
the isospin and the charge are decoupled completely, the contribution of isospin to the 
free energy vanishes. 

Weak coupling limit: In order to discuss the possibility of the existence of BEC, 
we consider the problem in weak coupling limit u —>■ 0. And isospin-isospin reaches its 
maximal correlation. At low temperature, however, we do not take string hypothesis for 
simplicity. Because limc^o -^n(a^) = together with eqs. (H) and eqs. (H), we obtain 

1 (3e"" - l)(e-"« + 1) 



pik) 



(25) 



2n (3e2^o + 1)(1 - e-=o) 
where Eq = (k'^—p) /T. The positive definition of p{k) require that the chemical potential 
is negative. As we known the density of scalar boson is 27rp = 1/ (1— e"^'') which prevents 
the BEC in ID and 2D system because of the infrared divergence. However, the density 
function ( p5D still does not resolve this problem. Consequently, BEC does not happen 
in this model. 

To summarize, we discussed the general thermodynamics of one dimensional SU(2) 
bosons with 5-function interaction by using the strategy of TBA. It was shown that 
the ground state is an isospin "ferromagnetic" state which differs from the ground 
state of ID fermions, while at high temperature, it is "paramagnetic" state and the 
magnetic property is dominated by Curie's law. In strong coupling limit, we obtain 
the exact expression of the dependence of chemical potential, entropy and specific heat 
on temperature which are Fermi-liquid like, while in weak coupling limit, we found the 
infrared divergence of charge roots density function prevents the existence of BEC. 

This work is supported by Trans-Century Training Program Foundation for the 
Talents and EYF98 of China Ministry of Education. SJG thanks D.Yang for helpful 
discussions. 
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